Classifying the Isolated Zeros of Asymptotic Gravitational Radiation by Tendex and 

Vortex Lines 
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A new method to visualize the curvature of spacetime was recently proposed. This method finds 
the eigenvectors of the electric and magnetic components of the Weyl tensor and, in analogy to 
the field lines of electromagnetism, uses the eigenvectors' integral curves to illustrate the spacetime 
curvature. Here we use this approach, along with well-known topological properties of fields on closed 
surfaces, to show that an arbitrary, radiating, asymptotically flat spacetime must have points near 
null infinity where the gravitational radiation vanishes. At the zeros of the gravitational radiation, 
the field of integral curves develops singular features analogous to the critical points of a vector 
field. We can, therefore, apply the topological classification of singular points of unoriented lines as 
a method to describe the radiation field. We provide examples of the structure of these points using 
linearized gravity and discuss an application to the extreme-kick black-hole-binary merger. 

PACS numbers: 04.30.-w, 04.20.-q, 02.40.Pc 



I. INTRODUCTION 

A recent study [l| proposed a method for visualiz- 
ing spacetime curvature that is well-suited for study- 
ing spacetimes evolved from initial data using numerical- 
relativity codes. The method first projects the Riemann 
curvature tensor R^ V(T p into a spatial slice, thereby split- 
ting it into two symmetric, trace-free spatial tensors, £ 
and 3 (see e.g. |2| and the references therein). These 
tensors are the spacetime-curvaturc analogs of the elec- 
tric and magnetic fields in Maxwell's theory. The electric 
tensor £ is familiar; it is the tidal field in the Newto- 
nian limit. The frame-drag field 3 (the magnetic cur- 
vature tensor) describes the differential frame dragging 
of spacetime. The eigenvectors of the tidal field provide 
the preferred directions of strain at a point in spacetime, 
and its eigenvalues give the magnitude of the strain along 
those axes. Similarly, the eigenvectors of the frame-drag 
field give preferred directions of differential precession of 
gyroscopes, and their eigenvalues give the magnitude of 
this precession [H, @, |j] ■ 

The study [l[ then proposed using the integral curves 
of these eigenvectors as a way to visualize the curvature 
of spacetime. Three orthogonal curves associated with 
£ , called tendex lines, pass through each point in space- 
time. Along each tendex line there is a corresponding 
eigenvalue, which is called the tendicity of the line. For 
the tensor 3, there is a second set of three orthogonal 
curves, the vortex lines, and their corresponding eigen- 
values, the vorticities. These six curves are analogous to 
the field lines of electromagnetism, and the six eigenval- 
ues to the electric and magnetic field strengths. The ten- 
dex and vortex lines, with their corresponding vorticities 
and tendicities, represent very different physical phenom- 
ena from field lines of electromagnetism; they allow one 
to visualize the aspects of spacetime curvature associated 
with tidal stretching and differential frame-dragging. In 
addition, each set of curves satisfies the constraint that 



its eigenvalues sum to zero at every point, since £ and 3 
are trace-free. 

Wherever the eigenvector fields are well-behaved, the 
tendex and vortex lines form extended, continuous fields 
of lines in a spatial slice. At points where two (or more) 
eigenvectors have the same eigenvalue, the eigenvectors 
are said to be degenerate. Any linear combination of the 
degenerate eigenvectors at these points is still an eigen- 
vector with the same eigenvalue; therefore, the span of 
these eigenvectors forms a degenerate subspacc. Singular 
features can appear at points of degeneracy, where many 
lines intersect, terminate, or turn discontinuously. The 
topology of unoriented fields of lines and their singular- 
points has been studied both in the context of general 
relativity and elsewhere. For example, Dclmarcelle and 
Hesselink Q studied the theory of these systems and ap- 
plied them to real, symmetric two-dimensional tensors. 
In the context of relativity, Penrose and Rindler Q ex- 
amined the topology of unoriented lines, or ridge systems, 
to characterize the principle null directions about single 
points in spacetime. Finally, Penrose p} also applied the 
study of ridge systems to human handprint and finger- 
print patterns. 

In this paper, we focus on the vortex and tendex lines 
and their singular points far from an isolated, radiating 
source. In Sec. HH we show that two of the vortex and 
tendex lines lie on a sphere (the third, therefore, is nor- 
mal to the sphere), and that the vortex and tendex lines 
have the same eigenvalues. Moreover, the two eigenval- 
ues on the sphere have opposite sign, and the eigenvalue 
normal to the sphere has zero eigenvalue. This implies 
that the only singular points in the lines occur when all 
eigenvalues vanish (i.e. when the curvature is exactly zero 
at the point, and all three eigenvectors are degenerate). 

In Sec. IIIII we employ a version of the Poincare-Hopf 
theorem for fields of integral curves to argue that there 
must be singular points where the curvature vanishes. 
Penrose, in a 1965 paper [8j, made a similar observation. 
There, he notes in passing that gravitational radiation 



must vanish for topological reasons, although he docs 
not discuss the point any further. Here we show that 
the topological classification of singular points of ridge 
systems can be applied to the tendex and vortex lines of 
gravitational radiation. This allows us to make a topo- 
logical classification of the zeros of the radiation field. 

In Sec. IIV1 we visualize the tendex and vortex lines 
of radiating systems in linearized gravity. We begin 
with radiation from a rotating mass-quadrupole mo- 
ment, the dominant mode in most astrophysical gravita- 
tional radiation. We then move to an idealized model of 
the "extreme-kick" configuration (an equal-mass binary- 
black-hole merger with spins antialigned in the orbital 
plane [9j ) . As we vary the magnitude of the spins in the 
extreme-kick configuration, we can relate the positions of 
the singular points of the tendex and vortex patterns to 
the degree of beaming of gravitational waves. We also 
visualize the radiation fields of individual higher-order 
multipole moments, which serve, primarily, as examples 
of patterns with a large number of singularities. Astro- 
physically, these higher multipoles would always be ac- 
companied by a dominant quadrupolc moment; we also, 
therefore, look at a superposition of multipoles. Since 
the tendex lines depend nonlincarly upon the multipoles, 
it is not apparent, a priori, how greatly small higher mul- 
tipoles will change the leading-order quadrupole pattern. 
Nevertheless, we see that for an equal-mass black- hole 
binary, higher multipoles make only small changes to the 
tendex line patterns. Finally, we discuss our results in 
Sec. H 

Throughout this paper we use Greek letters for space- 
time coordinates in a coordinate basis and Latin letters 
from the beginning of the alphabet for spatial indices 
in an orthonormal basis. We use a spacetime signature 
( — f- H — h) and a corresponding normalization condition 
for our tetrad. We will use geometric units, in which 
G = c=l, 

We will also specialize to vacuum spacctimes, where 
the Riemann tensor is equal to the Weyl tensor C^ up<J . 
To specify our slicing and to compute £ and B, we 
use a hypersurface-orthogonal, timelike unit vector, eo, 
which we choose to be part of an orthonormal tetrad, 
(eo, ei, e2, e^). We then perform a 3 + I split of the 
Weyl tensor by projecting it and its Hodge dual *C tlup<J = 






n.v a C a p ap into this basis, 



Sab — ^ajibv&Q eo , 

B a b = — *C ap b v ea IJ 'eo 1 ' . 



(1) 

(2) 



Here our convention for the alternating tensor is that 
£oi23 = +1 in an orthonormal basis. Note that, while the 
sign convention on B is not standard (see e.g. |H|), it has 
the advantage that £ and B obey constraints and evo- 
lution equations under the 3 + 1 split of spacetime that 
are directly analogous to Maxwell's equations in electro- 
magnetism 0, 3 . After the projection, we will solve the 
eigenvalue problem for the tensors £ and B in the or- 
thonormal basis, 



and we will then find their streamlines in a coordinate 
basis via the differential equation relating a curve to its 
tangent vector, 

ds 
Here s is a parameter along the streamlines. 



v"e a <- 



(4) 



II. GRAVITATIONAL WAVES NEAR NULL 
INFINITY 

Consider a vacuum, asymptotically flat spacetime that 
contains gravitational radiation from an isolated source. 
We are specifically interested in the transverse modes of 
radiation on a large sphere S near future null infinity. To 
describe these gravitational waves, we use an orthonor- 
mal tetrad (e ,e 1 ,e2,e 3 ), with e timelike and e2,e3 
tangent to the sphere, and we associate with this tetrad 
a corresponding complex null tetrad, 



1 



1 = — =(e + ei) 



« = -|( e o-ei), 



m = — = (e 2 + ie 3 ), m = — =(e 2 - ie 3 ) . (5) 

V2 V2 

Here, I is tangent to outgoing null rays that pass through 
S and strike a sphere at null infinity. We enforce that the 
null tetrad is parallelly propagated along these rays, and 
that it is normalized such that l^n^ = —m p fh^ = —1 
(all other inner products of the null tetrad vanish) . With 
these ra ys, we can associate Bondi-type coordinates (see 
e.g. pj],[l2|]) on a sphere at future null infinity with those 
on S. The timelike vector eo specifies our spatial slicing 
in this asymptotic region. When the orthonormal and 
null tetrads are chosen as in E q. (El ), £ and B are related 
to the complex Weyl scalars [13| . With the Newman- 
Penrose conventions appropriate to our metric signature 
(see, e.g., [10|). and our convention in Eq. ^f, one can 
show that 



'2*0 



Sab + iB a b = 



*3- 



*1 

- * 5 



«(*4 



-*s) 

-*o) 



^0 + 



>l"l 



(6) 



-4-2 



where * indicates entries that can be inferred from the 
symmetry of £ and B. 

In an asymptotically flat spacetime, the peeling the- 
orem [13| ensures that "J 4 ~ r _1 (with r an affine pa- 
rameter along the rays), and that the remaining Weyl 
scalars fall off with progressively higher powers of r, 
*3 ~ r~ 2 ,^ 2 ~ r~ 3 ,^i — r -4 , and * ~ ?*~ 5 • Asymp- 
totically, only $4 contributes to £ and B, 



Sab + iBab 



'0 
* 4 i^4 
.0 i* 4 -*4, 



(7) 



£ ab V = Xv a , 



(3) 



We see immediately that one eigenvector of both £ and B 
is the "radial" basis vector ei, with vanishing eigenvalue. 



The remaining 2x2 block is transverse and traceless, and 
the eigenvectors in this subspace have a simple analytical 
solution. The eigenvalues arc A± = ±|^4|/2 for both 
tensors, and the eigenvectors of £ have the explicit form 



v± 



-gggeg + (£22 - A±)e 3 

Im^ 4 e 2 + (Rc^ 4 =p |^ 4 |)e 3 
^(Im^^lRc^Tl^l) 2 



(8) 



The eigenvectors of 3 are locally rotated by ±7r/4 with 
respect to those of £ Q. As a result, although the global 
geometric pattern of vortex and tendex lines may differ, 
their local pattern and their topological properties on S 
will be identical. Moreover, when the eigenvalues of £ 
(the tendicity of the corresponding tendex line) vanish, 
so must those of 3 (the vorticity of the vortex lines). In 
arguing that the radiation must vanish, we can, therefore, 
focus on the tendex lines on S without loss of generality. 
Physically, however, both the vortex and the tendex lines 
are of interest. Similarly, since the two sets of tendex lines 
on S have equal and opposite eigenvalue and are orthogo- 
nal, we need only consider the properties of a single field 
of unoriented lines on S in order to describe the topo- 
logical properties of all four tendex and vortex lines on 
the sphere. Note that thus far we leave the coordinates 
(x 2 ,x 3 ) on S unspecified. We will assume that these 
coordinates are everywhere nonsingular, for instance by 
being constructed from two smooth, overlapping charts 
onS*. 



III. THE TOPOLOGY OF TENDEX PATTERNS 
NEAR NULL INFINITY 

Before investigating the properties of the tendex lines 
on S, we first recall a few related properties of vector 
fields on a 2-sphere. A well-known result regarding vec- 
tor fields on a sphere is the "hairy-ball theorem." This 
result states, colloquially, that if a sphere is covered with 
hairs at each point, the hair cannot be combed down 
everywhere without producing cowlicks or bald spots. 
The hairy-ball theorem is a specific illustration of the 
Poincare-Hopf theorem, applied to a 2-sphere. On a 2- 
sphere, this theorem states that the sum of the indices 
of the zeros of a vector field must equal the Euler char- 
acteristic, x, of the sphere, specifically \ = 2. The index 
of a zero of a vector field (also called a singular point) 
can be found intuitively by drawing a small circle around 
the point and traveling once around the circle counter- 
clockwise. The number of times the local vector field 
rotates counterclockwise through an angle of 27r during 
this transit is the index. More precisely, we can form a 
map from the points near a zero of the vector field to the 
unit circle. To do this consider a closed, oriented curve 
in a neighborhood of the zero, and map each point on 
the curve to the unit circle by associating the direction 
of the vector field at that point with a particular point on 





FIG. 1: Illustrations of the two types of half-index singulari- 
ties for ridge systems on a two-dimensional space. On the left 
is a loop singularity with index i = 1/2, and on the right is a 
triradius with i — —1/2. 



the unit circle. The index is the degree of the map (the 
number of times the map covers the circle in a positive 
sense). For the zero of a vector field, the index is a pos- 
itive or negative integer, because we must return to the 
starting point on the unit circle as we finish our circuit 
of the curve around the zero. 

The concept of an index and the formal statement of 
the Poincare-Hopf theorem generalizes naturally to ridge 
systems, fields of unoriented lines such as the tendex lines 
on S. For ridge systems on the sphere, the index of 
a singular point can be a half- integer J5|. Intuitively, 
this can occur because fields of lines do not have orien- 
tation. As one traverses counterclockwise about a small 
circle around a singular point, the local pattern of lines 
can rotate through an angle of ±7r during the transit. 
We illustrate the two fundamental types of singularity 
in Fig. [TJ which, following |||, we call loops for index 
i = 1/2 and triradii for i = —1/2. One can argue that 
the Poincare-Hopf Theorem holds for ridge systems, by 
noting that we can create a singular point with integer 
index by bringing two half-index singularities together 
(see Fig. [5] for a schematic of the creation of a singularity 
of index i = 1 from two loop singularities). Ridge pat- 
terns near singularities with integer index i = ±1 can be 
assigned orientations consistently; they must, therefore, 
have the same topological properties as streamlines of 
vector fields (which, in turn, have the same properties as 
the underlying vector fields themselves). By arguing that 
one can always deform a ridge system so that its singular 
points have integer index, one can see that the sum of 
the indices of a ridge system on a sphere must equal the 
Euler characteristic of the surface, % = 2 (see and the 
references therein for a more formal statement and proof 
of this theorem) . In Fig. [3] we show several other ridge 
singularities with integer index for completeness. In the 
top row, we show three patterns with index i = 1, and 
in the bottom left, we sketch a saddle type singularity 
with index i — ~ 1 . All of these patterns can be consis- 
tently assigned an orientation and, thus, have the same 
topological properties as vector field singularities. 

Having arrived at the result that the tendex lines on 
S must have singular points in a general, asymptotically 
flat vacuum spacetime, we now recall the fact that the 



singular points appear where there is a degenerate eigen- 
value of the tidal tensor. From the result of Sec. [D] 
the only degeneracies occur where the curvature van- 
ishes completely, and it follows therefore that there must 
be points of vanishing curvature on S. In general we 
would expect the radiation to vanish at a minimum of 
four points, as Penrose [8| had previously noted. In this 
case there would be four loop singularities with index 
i = 1/2, whose index sums to \ — 2. As we highlight 
in Sec. HVi where we show several examples of multipo- 
lar radiation in linearized theory, the number of singular 
points, the types of singularities, and the pattern of the 
tendex lines contains additional information. 

Additional symmetry, however, can modify the struc- 
ture of the singular points, as we see in the simple exam- 
ple of an axisymmctric, head-on collision of two nonspin- 
ning black holes. Axisymmctry guarantees that the Weyl 
scalar $4 is purely real when we construct our tetrad, Eq. 
([5]); by choosing e2 and e3 to be the orthonormal basis 
vectors of spherical polar coordinates on S, eg and e^ 
[14} . Using the relation \&4 = — h+ + ih x , we see that 
the waves are purely + polarized. By substituting this 
relationship into 0, we also see that eg and e^ are the 
eigenvectors whose integral curves are the tendex lines. 
The tendex lines, therefore, are the lines of constant lati- 
tude and longitude, and the singular points reside on the 
north and south poles of S. Their index must be i = 1, 
and the local pattern at the singularity will resemble the 
pattern at the top left of Fig. [3] for one set of lines, and 
the image on the top right of Fig. [3] for the other set (see 
also [J]). In this special situation, axisymmetry demands 
that there be two singular points on the axis, rather than 
four (or more). Moreover, these singular points are each 
generated from the coincidence of two loop singularities, 
with one singular point at each end of the axis of sym- 
metry. Similarly, if ^4 were purely imaginary, then the 
radiation would only contain the x polarization. The 
7r/4 rotations of the unit spherical vectors would then be 
the eigenvectors of the tidal field, and the two singular- 
ities at the poles would resemble that illustrated at the 
top middle of Fig. |3J 

It is even conceivable that four loops could merge 
into one singular point. This singularity would have the 
dipolelike pattern illustrated at the bottom right of Fig. 
[3j and it would have index i = 2. Though this situation 
seems very special, we show in the next section how a 
finely tuned linear combination of mass and current mul- 
tipolcs can give rise to this pattern. Because there is 
only one zero, the radiation is beamed in the direction 
opposite the zero, resulting in a net flux of momentum 
opposite the lone singular point. 

Before concluding this section, we address two possible 
concerns. The scalar ^4 = C^ vpcr n^fh v n p fh c ' depends 
both on the curvature and on the chosen tetrad. We 
first emphasize that the singular points we have discussed 
have nothing to do with tetrad considerations, in particu- 
lar with the behavior of the vectors tangent to the sphere, 
m, and fh. Though these vectors will also become singu- 




FIG. 2: An illustration of the formation of a singularity with 
index i = 1 from two loop singularities with index i = 1/2. 
The local structure of the two loops is shown in the top left, 
and the arrow represents, schematically, how they might join 
together into the extended pattern at the top right. Finally, 
the two loop singularities can be brought together until they 
coincide (which we represent by an arrow pointing to the im- 
age at the bottom). This resulting local pattern can be as- 
signed an orientation and is equivalent to the singular point 
of a vector field. 




FIG. 3: Top row: Diagrams of three orientable ridge patterns, 
which can be made from a combination of two loops, all with 
index i — 1. Bottom row: On the left is an orientable ridge 
pattern with index i — — 1 (which is identical to a saddle point 
of a vector field). It can be constructed by joining two triradii 
singularities. The figure on the right shows a dipolelike pat- 
tern with index i — 2, which can come from the coincidence 
of four loops. 



lar at points on the sphere, we are free to use a different 
tetrad on S in these regions, just as we can cover the 
sphere everywhere with smooth coordinates using over- 
lapping charts. Secondly, the vanishing of radiation does 
not occur due to the null vector n coinciding with a prin- 
ciple null direction of the spacetime. We note that, if ^4 
vanishes at a point on S, then a change of basis cannot 
make ^4 (or any of the other curvature scalars) nonva- 
nishing. For example, a rotation about I by a complex 
parameter a induces a transformation on the other basis 



vectors, 



I' = I, 

m! = m + al . 

fh' = fh + al , 

n' = n + am + am + aal . 



(9) 



Under this rotation, ^4 transforms as 

% = *4 + 4a* 3 + 6a 2 *2 +4a 3 *i +a 4 * , (10) 

which vanishes when the Weyl scalars are zero in the 
original basis. The remaining scalars transform analo- 
gously, and the other independent tetrad transformations 
are also homogeneous in the Weyl scalars (see e.g. |10j). 

IV. EXAMPLES FROM LINEARIZED GRAVITY 

We now give several examples of the tendex and vor- 
tex patterns on S from weak-field, multipolar sources. 
We first investigate quadrupolar radiation, produced by 
a time- varying quadrupolc moment. For many astro- 
physical sources, such as the inspiral of comparable mass 
compact objects, the gravitational radiation is predom- 
inantly quadrupolar. As a result, our calculations will 
capture features of the radiation coming from these as- 
trophysical systems. We will then study a combination 
of rotating mass- and current-quadrupole moments that 
are phase-locked. The locking of these moments was ob- 
served by Schnittman ct al. in their multipolar analysis 
of the extreme- kick merger [15[ . We conclude this sec- 
tion by discussing isolated higher multipoles. Although 
it is unlikely that astrophysical sources will contain only 
higher multipoles, it is of interest to see what kinds of 
tendex patterns occur. More importantly, while the tidal 
tensor is a linear combination of multipoles, the tendex 
lines will depend nonlincarly on the different moments. 
Actual astrophysical sources will contain a superposition 
of multipoles, and it is important to see how superposi- 
tions of multipoles change the leading-order quadrupolc 
pattern. 

We perform our calculations in linearized theory about 
flat space, and we use spherical polar coordinates and 
their corresponding unit vectors for our basis. One can 
compute from the multipolar metric in jig ] that for a 
symmetric, trace-free (STF) quadrupolc moment I a b, the 
leading-order contributions to £ and on S are 



B 



e=2) 



1 

2~7 

1 

r 



(4) Z afc (t -r) + £ac (4) lcd(t ~ r)e db 



TT 



c(a ( % c (t-r) 



TT 



(11) 

(12) 



Here, the superscript ( 4 ) indicates four time derivatives, 
TT means to take the transverse-traceless projection of 
the expression, and e ac is the antisymmetric tensor on 



a sphere. In this expression, and in what follows, the 
Latin indices run only over the basis vectors eg and e^, 
and repeated Latin indices are summed over even when 
they are both lowered. 



A. Rotating Mass Quadrupole 

As our first example, we calculate the STF quadrupole 
moment of two equal point masses (with mass Af/2) sepa- 
rated by a distance a in the equatorial plane, and rotating 
at an orbital frequency il. Wc find that 

w l ge {t-r) = Ma 2 fl 4 (l + cos 2 9) cos{2[0 - Q(t - r)]} , 
{i) l g(t> (t - r) = -2Ma 2 n i cos6sm{2[<f)-n(t-r)]}, 
%(t-r) = -^l ee {t-r). (13) 

By substituting these expressions into Eqs. (TTT1) and (J5J, 
we find the eigenvectors of the tidal field. We can then 
calculate the tendex lines on the sphere by solving Eq. ^ 
with a convenient normalization of the parameter along 
the curves, 



d0 1 (4)T 

as r 



1 



(Is 



-(-^lee - A+) . 



(14) 
(15) 



Here, A+ is the positive eigenvalue. The differential equa- 
tion for the vortex lines [found from the corresponding 
frame-drag field of Eq. (TT2|) ]. has the same form as those 
of the tendex lines above; however, one must replace 
^Igrp in the first equation by ^Igg and ^Igo by —^Ig^ 
in the second equation. 

We show the tendex and vortex lines corresponding to 
the positive eigenvalues in Figs. |4] and respectively, at 
a retarded time t — r = 0. We also plot the magnitude 
of the eigenvalue on the sphere, using a color scheme in 
which purple (darker) regions at the poles correspond 
to large eigenvalues and yellow (lighter) colors near the 
equator are closer to zero. Both the tendex and vortex 
lines have four equally spaced loop singularities on the 
equator at the points where the field is zero (the two on 
the back side of the sphere are not shown). Because the 
vortex and tendex lines must cross each other at an angle 
of 7r/4, the global geometric patterns are quite different. 

We note here that these two figures also provide a visu- 
alization for the transverse-traceless, "pure-spin" tensor 
spherical harmonics [16| . For example, we can see that 
the mass-quadrupolc tendex lines are the integral curves 
of the eigenvectors of the real part of the £ = 2, m = 2 
electric- type transverse-traceless tensor harmonic. First, 
the tendex lines correspond to the electric- type harmonic, 
because the tidal tensor is even under parity. Second, the 
radiation pattern will not contain an £ = 2, m = har- 
monic, because the overall magnitude of the quadrupolc 
moment of the source is not changing in time; also, the 
£ = 2, m = ±1 harmonics arc absent because the source 



is an equal-mass binary and is symmetric under a rota- 
tion of 7r. Finally, the I = 2, m = — 2 moment is equal 
in magnitude to the m = 2 harmonic, since the tidal ten- 
sor is real. By similar considerations, we can identify the 
vortex lines of the mass quadrupole as a visualization of 
the real parts of the 1 = 2, m = 2 magnetic- type tensor 
harmonics. 

In addition, the eigenvalue (the identical color patterns 
of Figs. [4] and [5]) is given by the magnitude of the sum of 
spin- weighted spherical harmonics, 



A, 



\-2Y0: 



-•211-<i 



(16) 



One can see this most easily by using the symmetries 
described above, the expression for the eigenvalue A+ = 
I SI/4 1 /2, and the spin-weighted spherical harmonic decom- 
position of ^4. It is also possible to verify this expres- 
sion using the tensor harmonics above and the standard 
relations between tensor spherical harmonics and spin- 
weighted spherical harmonics (see e.g. [III). Radiation 
from numerical spacetimes is usually decomposed into 
spin- weighted spherical harmonics, and, as a result, the 
pattern of the eigenvalue is familiar. The tendex lines, 
however, also show the polarization pattern of the waves 
on S (a feature that numerical simulations rarely explic- 
itly highlight). FigureS](and the accompanying negativc- 
tendicity lines not shown) gives the directions of preferred 
strain on S, and hence the wave polarization that can 
be inferred from gravitational-wave-interferometer net- 
works such as LIGO/VIRGO. Thus, visualizations such 
as Fig. 0] give complete information about the gravita- 
tional waves passing through S. 



Rotating Mass and Current Quadrupoles in 
Phase 



As our second example, we will consider a source that 
also has a time- varying current-quadrupole moment, S a b- 
In linearized theory, one can show that the tidal tensor 
and frame-drag field of a current quadrupole are simply 
related to those of a mass quadrupole. In fact, B a b of 
the current quadrupole has exactly the same form as £ a b 
of a mass quadrupole, Eq. (fTTj) . when one replaces ^X a b 
by (4/3p 4 )«S a &- Similarly, £ a b of the current quadrupole 
is identical to B a b of a mass quadrupole, Eq. (fT2"j) . when 
^l ab is replaced by -(4/3)( 4 )<S a6 . 

We impose that the source's mass- and current- 
quadrupole moments rotate in phase, with frequency ft, 
and with the current quadrupole lagging in phase by n/2. 
This arrangement of multipoles models the lowest mul- 
tipolcs during the merger and ringdown of the extreme- 
kick configuration (a collision of equal- mass black holes in 
a quasicircular orbit that have spins of equal magnitude 
lying in the orbital plane, but pointing in opposite direc- 
tions), when the mass- and current-quadrupole moments 
rotate in phase [15| • The relative amplitude of the mass- 
and current-multipoles depends upon, among other vari- 
ables, the amplitude of the black-holes' spin. We, thcrc- 




FIG. 4: (color online). The positive tendex lines on S gener- 
ated by a rotating quadrupole moment in linearized gravity. 
The sphere is colored by the magnitude of the eigenvalue; pur- 
ple (darker) areas at the poles corresponding to a large eigen- 
value, and yellow (lighter) areas near the equator indicate a 
value closer to zero. Four loop singularities appear equally 
spaced on the equator at the points of vanishing tendicity. 



fore, include a free parameter C in the strength of the 
current quadrupole which represents the effect of chang- 
ing the spin. An order-of-magnitude estimate based on 
two fast-spinning holes orbiting near the end of their in- 
spiral indicates that that their amplitudes could be nearly 
equal, C — 0(1). To determine the exact relative am- 
plitude of the mass- and current-quadrupole moments of 
the radiation would require comparison with numerical 
relativity results. 

We calculate the current-quadrupole moment by scal- 
ing the mass quadrupole by the appropriate factor of C 
and letting the term 2[<f> — fl(t — r)] in the equations for 
Z a b(t — r) become 2[<j>— fl(t — r)] —tt/2 in the correspond- 
ing expressions for S a b(t — r). In linearized theory, the 
tidal tensor and frame-drag fields of the different mul- 
tipoles add directly. As a result, the equations for the 
tendex lines have the same form as Eqs. (fl4"|) and (|15|) . 
but one must now replace the mass quadrupole ^'Xg^ by 
( 4 >Terp — (4/3) ( 4 ) See in the first expression and ^Igg by 
^l eg + (4/3) ^S^ in the second. 

First, we allow the current quadrupole to be half as 
large as the mass quadrupole, C = 1/2. We show the 
positive tendex lines and positive eigenvalue in Fig. [5] 
Because of the relative phase and amplitude of the two 
moments, the tensors add constructively in the northern 
hemisphere and destructively in the southern hemisphere 
on S. This is evident in the eigenvalue on the sphere in 
Fig. [5J which, one can argue, is now given by an unequal 





FIG. 5: (color online). As in Fig. SI we show the positive 
vortex lines and their magnitude of the eigenvalue on S (using 
the same coloring as in that figure) . The loop singularities lie 
at the same locations as they do for the tendex lines, but they 
are locally rotated by 7r/4. 



superposition of spin-weighted spherical harmonics, 
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with b < 1. As in previous figures, dark colors (black 
and purple) represent where the eigenvalue is large, and 
light colors (white and yellow) show where it is nearly 
zero. While the singular points are still equally spaced 
on a line of constant latitude, they no longer reside on 
the equator; they now fall in the southern hemisphere. 
This is a direct consequence of the beaming of radiation 
toward the northern pole. 

The case shown above has strong beaming, but it is 
possible to make the beaming more pronounced. To 
get the greatest interference of the multipoles, the mass 
and current quadrupoles must have equal amplitude in 
the tidal field. Because the tidal field of the current 
quadrupole is 4/3 as large as the tidal field of the mass 
quadrupolc, setting C — 3/4 gives the strongest con- 
structive interference in the tidal fields. In this case, the 
eigenvalue vanishes at just one point, the south pole, and 
the eigenvalue can be shown to be proportional to just a 
single spin- weighted spherical harmonic, 



A. 
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As a result, the four equally spaced singular points of the 
tendex lines must coincide at one singular point whose 
index must be i = 2. This is precisely the dipolelike 
pattern depicted in Fig. [3l We show the tendex lines 
around the south pole in Fig. [7J The vortex lines are 
identical to the tendex lines, but they are globally rotated 
by 7r/4 in this specific case. 



FIG. 6: (color online). The positive tendex lines on S gen- 
erated by the superposition of rotating mass- and current- 
quadrupole moments, n/2 out of phase, in linearized gravity. 
The sphere is colored by the tendicity as in Fig. [4] Interfer- 
ence between the moments leads to beaming of the radiation 
toward the north pole. Similarly, the singular points of the 
tendex lines now fall on a line of constant latitude in the 
southern hemisphere. 



We see that the beaming can be maximized by care- 
fully tuning the phase and amplitude of the mass- and 
current-quadrupole moments. Interestingly, the maxi- 
mally beamed configuration corresponds with the coin- 
cidence of all singular points at the south pole in the 
radiation zone. Whether this degree of beaming could 
occur from astrophysical sources is an open question. 



C. Higher Multipoles of Rotating Point Masses 

We also investigate the effect of including higher mul- 
tipoles on the tendex lines on S. For the orbiting, non- 
spinning, point masses of the first example the next two 
lowest multipoles arise from the current octopole (the 
I = 3 STF moment [ljl) and the mass hexadecapole (th e 
I = 4 STF moment). From the multipolar metric in [la |. 
one can show that the tidal field for these two moments 
are 
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where the index r indicates contraction with the radial 
basis vector e r , and so repeated r indices do not indi- 
cate summation. The STF current-octopolc moment, can 
be expressed compactly as Sijk — (L z N x :, A x\) STF , where 




^l^ rr (t-r) = - i6) lg grr {t - r) . 



(22) 



FIG. 7: (color online). South polar region of the tendex line 
pattern of a gravitational wave generated by rotating mass- 
and current-quadrupole moments. The amplitude and phase 
of the moments are chosen so that the radiation vanishes only 
at the south pole. The purple (darker) areas indicate a pos- 
itive eigenvalue, while yellow (lighter) areas are values closer 
to zero. The singularity at the pole has index i — 2. 



L % N is the Newtonian angular momentum and x A is the 
position of one of the point masses. The superscript 
STF indicates that all indices should be symmetrized, 
and all traces removed. In Cartesian coordinates, the 
vectors have the simple forms L^r = (0, 0, Mav/4) and 
x^ = (a/2)(cos[Qt],sin[Q£],0), where O is the Keple- 
rian frequency and v is the relative velocity. Similarly, 
one can write the STF mass-hexadecapolc moment as 
lijkl = M (x 1 A x 3 A x\x l A ) STF , for the same vector x A as 
above. Because these tensors have many components, 
we shall only list those that are relevant for finding the 
tendex lines. We will also define a = <fi — fi(£ — r) for 
convenience. For the current octopole the relevant com- 
ponents are 

{5) S 99r (t~r) = (5 cos 9 + 3 cos 30) sin 2a, 

Ma 3 vQ 5 



(c\ _ , lviavii 

Se<j>r(t — r) = cos2#cos2a. 

&S^ r (t-r) = ~^S 9 er(t-r), 
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and for the mass hexadecapole they are 

■ [(cos 2 9 + cos 49) cos 2a 
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v 'Xg^rr (t — r) = [cos W sin 2a 

- 128 sin 2 9 cos 9 sin 4a] , 



The tendex lines of the current octopole can be found by 
solving the system of differential equations in Eqs. (|14|) 
and (JTSJ) by substituting ^X e4> by (5) <W/2 and ^l e9 
by — ( 5 )<Se^ r /2 Similarly, for the mass hexadecapole, one 
must make the substitutions of ' 4 'le0 by ^>Ig^ rr /Yl and 
( A >lee by ( 6 >Igg rr /12 in the same equations. 

In Fig. [8] we show the tendex line pattern for the cur- 
rent octopole, and in Fig. [H] we show the pattern for the 
mass hexadecapole. Together with the mass quadrupolc, 
Fig. |4l these are the three lowest multipole moments 
for the equal-mass circular binary. For the current oc- 
topole, there are eight triradius singular points and 12 
loop singularities (and thus the net index is two). Four 
of the loop singularities remain equally spaced on the 
equator, at the same position of those of the quadrupole, 
but the remaining singularities appear at different points 
on S. The mass hexadecapole has eight loop singulari- 
ties equally spaced on the equator, and there are integer- 
index saddlc-point-like singularities at each pole. 

Gravitational radiation from astrophysical sources will 
likely not be dominated by these higher multipoles. Nev- 
ertheless, these figures are of interest as examples of ten- 
dex lines with many singular points and as visualizations 
of tensor harmonics. By analyzing the symmetries in a 
way analogous to that discussed in Sec. IIV A[ we can 
identify the currcnt-octopolc tendex lines with the inte- 
gral curves of £ = 3, m — 2 magnetic- type harmonics, 
and we can associate the mass-hexadecapole lines with 
those of the £ = 4, m = 4 electric-type harmonics. In 
the case of the mass hexadecapole, the m = 2 moment is 
not ruled out by symmetry, but it is suppressed relative 
to the m = 4 moment. This occurs because the m = 4 
moment oscillates at twice the frequency of the m = 2 
moment, and the tidal tensor for this higher-order mo- 
ment is given by taking six time derivatives of the STF 
moment, Eq. (f2"2")l . This enhances the m = 4 radiation 
by a factor of 2 6 over the m = 2 contribution. Sim- 
ilarly, we can relate the eigenvalue to the magnitude of 
the corresponding sum of s = — 2 spin- weighted spherical 
harmonics, and the tendex line patterns to the the polar- 
ization directions that could be inferred from networks 
of gravitational- wave interferometers. 

Finally, we show the pattern generated from the lin- 
ear combination of the three lowest multipole moments 
in Fig. [TO] Any astrophysical source will contain several 
multipoles, with the quadrupolc being the largest. The 
tendex lines depend nonlincarly on the multipoles, and 
it is important, therefore, to sec to what extent higher 
multipoles change the overall pattern. We find the to- 
tal tidal tensor by linearly combining the tidal tensor 
of each individual moment, and we then find the eigen- 
vectors and tendex lines of the total tidal tensor. The 
pattern formed from the combination of multipoles de- 
pends upon the parameters of the binary; in making this 
figure we assumed (in units in which M = 1) a separa- 
tion of a = 15, an orbital frequency O = a~ 3 ' 2 , and a 





FIG. 8: (color online). The tendex lines of a current-octopole 
moment of an equal-mass, circular binary of point masses. 
The colors on the sphere represent the tendicity, with the 
same scale described in Fig. [4] The current octopole also has 
four loop singularities on the equator (at the same position of 
those of the rotating quadrupole), but it has eight additional 
loops and eight triradius singularities off of the equator. Only 
half of the singular points are visible on the sphere; the other 
half appear on the back side. 



FIG. 9: (color online). The tendex lines on S of the mass hex- 
adecapole of an equal-mass, circular binary of point masses, 
with the sphere colored by the tendicity as in Fig. [4] The hex- 
adecapole has eight loop singularities equally spaced on the 
equator and two saddle-point-like singularities (from the coin- 
cidence of two triradius singularities at a point) at the poles. 
Again, only half are visible in the figure. Four of the singular 
points on the equator coincide with those of the quadrupole, 
but the remaining four appear at different locations. 



velocity v = afl = a -1 / 2 . When these higher moments 
are combined with the mass quadrupole, the tendex line 
structure resembles that of the mass quadrupole. The 
pattern is deformed slightly, however, by the presence 
of the higher multipoles. The loop singularities on the 
equator are no longer evenly spaced; rather, the pair il- 
lustrated (and the corresponding pair which is not visi- 
ble) are pushed slightly closer together. 



V. CONCLUSIONS 

Tendex and vortex lines provide a new tool with which 
to visualize and study the curvature of spacetime. Funda- 
mentally, they allow for the visualization of the Riemann 
tensor, through its decomposition into two simpler, trace- 
free and symmetric spatial tensors. These tensors, £ and 
6, can be completely characterized by their eigenvectors 
and corresponding eigenvalues. The integral curves of 
these eigenvector fields are easily visualized, and their 
meaning is well understood; physically, the lines can be 
interpreted in terms of local tidal strains and differential 
frame-dragging. Here, the simple nature of these lines 
allows us to apply well-known topological theorems to 
the study of radiation passing through a sphere near null 
infinity. 

Tendex line patterns must develop singularities (and 
thus have vanishing tendicity) on a closed surface. When 



we applied this fact to the tendex lines of gravitational 
radiation near null infinity from arbitrary physical sys- 
tems, we could easily show that the gravitational radi- 
ation must at least vanish in isolated directions. Al- 
though this result is somewhat obvious in retrospect and 
has been noted before Q, the result does not appear to 
be well-known. We also began exploring the manner in 
which these singular points can provide a sort of finger- 
print for radiating spacctimes. The essential elements 
of this fingerprint consist of the zeros of the curvature 
on the sphere, together with the index and the tendex 
line pattern around these zeros. We studied these pat- 
terns for a few specific examples, such as the four equally 
spaced loops of a rotating mass quadrupole. A more in- 
teresting case is that of a radiating spacetime composed 
of locked, rotating mass and current quadrupoles, which 
can be thought of as a simplified model of the late stages 
of the extreme-kick black-hole-binary merger. Here, the 
shifted positions of the singular points of the tendex pat- 
tern provide a direct illustration of gravitational beaming 
for this system. By seeking the most extreme topologi- 
cal arrangement of singular points, we also described a 
maximally beaming configuration of this system. 

The radiation generated by higher-order STF multi- 
pole moments gives more complex examples of tendex 
and vortex patterns, with many singular points of varied 
types. Additionally, we argued that their tendex and vor- 
tex patterns provide a visualization of the tensor spher- 
ical harmonics on the sphere; the eigenvalue illustrates 
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FIG. 10: (color online) . The tendex lines of a superposition of 
mass-quadrupole, current-octopole, and mass-hexadecapole 
moments of an equal-mass circular binary. It assumes a to- 
tal mass M — 1, a separation a = 15, an orbital frequency 
Q = a~ 3 ' 2 and a velocity v = afi = a -1 ' 2 . The sphere is 
colored by the tendicity in an identical way to that of Fig. 
IU When the tidal tensors of the three multipoles are com- 
bined, the net pattern is dominated by the quadrupole and 
contains only the four loops. The loop singularities are no 
longer equally spaced on the equator; the two pairs are pushed 
closer together due to the influence of the higher multipoles. 



the magnitude of these harmonics, and the lines show the 
tensor's polarization in an intuitive manner. The sum of 
the three multipoles illustrated in Fig. [TU] shows how in- 
cluding higher-order multipoles slightly deforms the pat- 
tern of quadrupole radiation to make a more accurate 



total radiation pattern of the equal-mass binary. Simi- 
lar illustrations of complete radiation patterns could be 
readily produced from numerical spacetimes, when ^4 is 
extracted asymptotically using a tetrad with appropriate 
peeling properties. Such visualizations, and their evolu- 
tion in time, could provide a useful method for visualizing 
the gravitational emission from these systems. 

This study of the tendex and vortex lines (and their 
singularpoints) of asymptotic radiation fields is one of 
several [4j exploring and developing this new perspective 
on spacetime visualization. Naturally, it would be of in- 
terest to extend the two-dimensional case here to a larger 
study of the singular points in the full, three-dimensional 
tendex and vortex fields. Methods to find and visual- 
ize the singular points (and singular lines) of 3D tensors 
have been discussed preliminarily in [17| . though there 
is still room for further work. We suspect that singular 
points will be important in visualizing and studying the 
properties of numerical spacetimes with these methods. 
Further, we expect that there is still much to be learned 
from the study of the vortexes and tendexes of dynamical 
spacetimes. 
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